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Abstract. We study pion condensation and the phase structure in a two- 
flavour Nambu-Jona-Lasinio model in the presence of baryon chemical potential 
fi and isospin chemical potential fij at zero and finite temperature. There is 
a competition between the chiral condensate and a Bosc-Einstcin condensate of 
charged pions. In the chiral limit, the chiral condensate vanishes for any finite 
value of the isospin chemical potential, while there is a charged pion condensate 
that depends on the chemical potentials and the temperature. At the physical 
point, the chiral condensate is always nonzero, while the charged pion condensate 
depends on fij and T. For T = /J, = 0, the critical isospin chemical potential 
/ij for the onset of Bose-Einstein condensation is always equal to the pion mass. 
For fi = 0, we compare our results with chiral perturbation theory, sigma-model 
calculations, and lattice simulations. 

Finally, we examine the effects of imposing electric charge neutrality and weak 
equilibrium on the phase structure of the model. In the chiral limit, there is a 
window of baryon chemical potential and temperature where the charged pions 
condense. At the physical point, the charged pions do not condense. 



PACS numbers: 14.40.Aq, 11.30.Qc, 12.39.-x, 21.65.-f 
1. Introduction 

There has been a tremendous effort in recent years to map out the phase diagram of 
QCD as a function of temperature and baryon chemical potential [H El 02 HI [5j [6l [7] . 
It is generally accepted that one can calculate the properties of strongly interacting 
matter at asymptotically high temperature or at asymptotically high densities using 
perturbative QCD. 

At sufficiently high density and low temperature, we know that QCD is in 
the colour-flavour locked (CFL) phase. This state is a colour superconducting 
state because the quarks form Cooper pairs in analogy to electrons in an ordinary 
superconductor. In this case, the original symmetry group of QCD, SU (3) c x SU (3)^ x 
SU(3)n x U(1)b, is broken down to SU(3) c +l+r which is a linear combination of the 
generators of the original group. This linear combination locks rotation in colour 
space to rotations in flavour space and this has given the name to the phase. In 
the CFL phase there is an octet of Goldstone modes which arises from the breaking 
of chiral symmetry and a singlet arising from the breaking of the baryon-number 
conserving group U(1)b- The CFL phase is a superfluid due to the breaking of the 
global U(1)b symmetry. This is analogous to the superfluidity encountered in in 
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condensed-matter systems such as 4 He. At very high densities, all nine modes are 
effectively massless since one can ignore the quark masses. This implies that the 
low-energy properties of the CFL phase can be described in terms of an effective 
field theory for the massless mesons [5J HI [TOl [TT] . At moderate densities, one cannot 
neglect the quark masses and chiral symmetry is broken explicitly. Thus only the 
superfluid mode is exactly massless, while the the other mesons acquire masses. 
The lightest massive modes are expected to be the charged and neutral kaons and 
if the chemical potentials are large enough, there is a transition to a Bose condensed 
phase. Bose condensation of kaons in the CFL phase has been studied in detail in 

Refs. Ha na us nsi nsi nzi nsi iia isa i5u . 

QCD at finite baryon chemical potential /is is not accessible by Monte Carlo 
simulations due to the complex fermion determinant. This is in contrast to QCD at 
finite isospin chemical potential \i\ (still at zero /is) where lattice simulations are 
possible since the functional determinant is real. Thus this is a system whose phase 
diagram one can study on the lattice as a function of a conserved charge. 

Lattice simulations [23J [24J [25] suggest that there is a deconfinement transition 
of pions at high temperature and low density, and Bose-Einstein condensation of 
charged pions at high isospin density and low temperature. In fact, the deconfinement 
transition and the transition to a charged pion condensate seem to coincide. 
The deconfinement transition is found by measuring the Polyakov loop and the 
measurements show a sharp increase (indicating deconfinement) at approximately the 
same temperature as the onset of pion condensation. 

Pion condensation and the phase diagram of two-flavour QCD have been in- 
vestigated using chiral perturbation theory [26l [2Zl EE], ladder QCD [29], the 
chiral quark model [3U], the linear sigma models [3TJ [3H [331 (Ml [351 NJL 
models [331 E3 EH ES SOI SD S21 S3, and Polyakov-loop NJL models [H [45]. The 
PNJL models suggest that deconfinement and the onset of Bose-Einstein condensation 
are two different transitions. Very recently, Abuki et al [46] have considered the 
possibility of probing the phase diagram of electrically neutral QCD at finite isospin 
chemical potential using an equilibrated gas of neutrinos. They found that a 
condensate of charged pions arises at large enough neutrino densities and small baryon 
densities, and that at even larger neutrino densities there is condensation of charged 
kaons as well. 

We note in passing the similarity between three-colour QCD at finite /i/ and 
two-colour QCD at finite [Ib E3j- The correspondence is given by identifying /i//2 
with [is-, the charged pion condensate with the diquark condensatey, and the isospin 
density with the quark density. 

Bose-condensed states or colour-superconducting states may be found in the 
interior of compact stars if the density is high enough. In contrast to hadronic matter 
in heavy-ion collisions, bulk matter in compact stars must (on average) be electrically 
neutral and so a neutrality constraint must be imposed [47, 45]. Similarly, bulk matter 
must be colour neutral and if the system is in a colour superconducting phase, one 
sometimes has to impose this constraint explicitly. It is automatically satisfied if one 
uses the QCD Lagrangian, but this is not so if one describes the system using NJL- 
type models. This is due to the fact that there are no gauge fields in this model and 
the SU(N C ) colour symmetry is global [15 ] W§ \ I5U 1 1ST] . 

f The diquark condensate in two-colour QCD does not break any local symmetries, only global ones. 
The system is therefore a superfluid but not a colour superconductor. 
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The advantage of models with quarks as microscopic degrees of freedom, such 
as the NJL model, is that one can investigate simultaneously the effects of finite 
baryon chemical potential and isospin chemical potential. In the present paper, we 
consider the two-flavour NJL model at finite baryon chemical potential and isospin 
chemical potential. We compute the phase diagram at zero and finite temperature as 
a function of these variables. We restrict ourselves to sufficiently low values of the 
baryon chemical potential such that there are no superconducting phases [48] [52] . 
We also investigate the effects on the phase diagram by imposing electric charge 
neutrality and /3-equilibrium. Our work is a generalization of the papers by Ebert 
and Klimenko 38, 39J to finite temperature and finite pion mass. 

The article is organized as follows. In Sec. [H we discuss the Lagrangian and 
the gap equations of NJL model. In Sec. [3j we discuss the phase diagram at zero as 
well as finite temperature. In Sec. IH we discuss the issues of charge neutrality and (3- 
equlibrium. In Sec.0 we investigate the phase diagram at zero and finite temperature 
imposing charge neutrality and /3-equlibrium. In Sec. [51 we summarize and conclude. 

2. Lagrangian and Gap equations 

In this section, we discuss the properties of the Lagrangian of the two-flavour NJL 



model. The Lagrangian can be written as [53J 

£ = Co + Ci + £ 2 , (1) 
where the various terms are 

£ = - m o)V- , (2) 

d = Gx [(<AV0 2 + (^) 2 + (^75</0 2 + (^75^) 2 ] , (3) 

£2 = G 2 [(^) 2 - (^rV>) 2 - (^ 75 V0 2 + (^75TV) 2 ] ■ (4) 



Here, m is the quark-mass matrix, which is diagonal in flavour space and contains 
the bare quark masses m u and mj. Moreover, r = (r 1; t%, r 3 ) where Tj (i = 1, 2, 3) are 
the Pauli matrices. G\ and G 2 are coupling constants. The quark field ip is an isospin 
doublet 




In the following, we take m u = ma- The Lagrangian (fTJ) has a global SU(N C ) symmetry 
as well as a U(1)b baryon symmetry. The latter reflects baryon number conservation. 
In the chiral limit, the Lagrangian (fTJ) has an SU(2)l x SU(2)r symmetry. Away from 
the chiral limit, this symmetry is reduced to SU(2)l + ji isospin symmetry. £\ has an 
additional U(1)a axial symmetry. £2 is 't Hooft's instanton-induced interaction term 
and breaks explicitly the U(1)a axial symmetry of £\ [54j[55] . In the following, we shall 
limit ourselves to study the standard NJL Lagrangian by choosing G\ — G2 = G/2 [53] 
and so Eq. (fTJ) reduces to 

£ = ^(ij^d^ - m )'/' + G [(ipip) 2 + (V^tVj) 2 ] ■ (6) 

We can characterize the system described by the Lagrangian ^ by the expectation 
values of the different conserved charges associated with the continuous symmetries. 
For each conserved charge Qi, we introduce a chemical potential fii. Note, however, 
that it is possible to specify the expectation values of different charges simultaneously 
only if they commute. In the present case, we introduce a chemical potential \ib 
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associated with the U{1)b baryon symmetry, as well as a chemical potential \ij 
associated with the third component of the SU(2)l + r isospin group. This is done 
by adding to the Lagrangian ©, the terms 

C B = ti B ^l°B4> , (7) 

Ci = nitfh*l> , (8) 

where B — diag(l/3, 1/3) and I3 = T$/2. We can then write the Lagrangian as 

C = j [ij^ - m + tij° + 6 m °T 3 ] i>+G [(^) 2 + (^ 7 5 ^) 2 ] , (9) 

where we have defined the quark chemical potential /i as well as Sfj, by 

M=f, (10) 

In the remainder of the paper we assume /j, > and <5/j > for simplicity. In terms of 
the chemical potentials for the u and the d-quarks, [i u and /id, the chemical potentials 
H and 8\x can be written as 

M = 2 G«« + /-*<*) , (12) 

Sf^ = ^(fhi- fid) • (13) 

In the chiral limit, the inclusion of the isospin chemical potential breaks the SU(2)l x 
5Z7(2)fl-symmetry of the Lagrangian to U(1)l x J7(1)h. At the physical point, it 
breaks the S , [/(2)i+fl-symmetry down to U(\)l+r- 

From the path-integral representation of the free energy density 



e 

the expression for the charge density Qi associated with the chemical potential ^ can 
be written as 

Qi=-JT- (15) 
We next introduce the auxiliary fields a and 7Tj by 

a = - 2G$ip , (16) 

7Tj = - 2Gtpij 5 Ti^p . (17) 
The Lagrangian §§§ can now compactly be written as 

C = ip [i^dft -m + (jsy + <5m7° t 3 - o - ij 5 TT a T a ] ip 

-^(^ 2 +^a) • (18) 

The original Lagrangian ([9]) can be recovered by using the equations of motion for 
the auxiliary fields a and 7Tj to eliminate them from Eq. (|18[) . The Lagrangian is now 
bilinear in the quark fields and so we can integrate them out exactly. We then obtain 
the following effective action for the composite fields a and ~K{ 

Sc« = ~ -N c Tr log [i'fd^ -m + M7° + Sfx-y°r 3 -a- i-/ 5 n a T a ] 



~~ I Jo (°* + na7Ta ^ ' (19) 
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where Tr denotes the trace and is over Dirac indices as well as space-time. 

We next introduce a nonzero expectation value for the fields a and 7Ti to allow 
for a chiral condensate and a charged pion condensate. The fields are then written as 



where a and tt\ are quantum fluctuating fields. In the mean-field approximation, we 
neglect the fluctuations of the quantum fields a and 7ri in the functional determinant. 
This approximation coincides with the leading order of the l/-/V c -expansion, where N c 
is the number of colours @. 

For notational simplicity, we introduce the quantities M and p given by 



where M is the constituent quark mass. Note that in the chiral limit, the chiral 
condensate breaks the SU(2)l x SU(2)r symmetry spontaneously down to SU(2)l + h 
in the usual manner. Moreover, we can always use the remaining [/(l)-symmetry to 
rotate away any nonzero value of (')/>7 E Tziji) ■ 

Note that we have introduced a single chiral condensate and not separate chiral 
condensates (uu) and {dd) for the u and the d quarks. We have chosen G\ = Gi and 
in this case the effective action (TT51) only depends on the sum of these condensates. 
In fact, it is easy to show that the two chiral condensates must be equal. This is in 
contrast to the calculations of Barducci et al j37], where the authors choose G\ — G/2 
and G2 = 0. In that case the effective action is not symmetric under permutation of 
the two chiral condensates and they are different. Our choice is motivated by the fact 
that U(1)a symmetry is broken in QCD due to instanton effects [5B]. In fact, setting 
G\ — G2 ^ in the present model means that axial symmetry is maximally violated. 
When G2 — 0, on the other hand, the Lagrangian Jl} is {/(l^-symmetric. 

Also note that we take the condensates M and p to be spacetime independent. 
Though it has been demonstrated that the phase diagram may contain phases where 
the condensates are non- uniform [57], also within the framework of the two-flavour 
NJL model [551 USB EES EH > taking this possibility into account is beyond the scope 
of this study. Still, it would be interesting to see how these non-uniform phases are 
affected by imposing neutrality constraints. 

Using standard techniques to evaluate the trace and using ft = —S e fi/0V, where 
is the thermodynamic potential and V is the volume of the system, we obtain 



7T! = - 2G{ipi^T 1 tp) + TTl , 



(20) 
(21) 



M = m - 2G(^ip) , 
p = - 2G(& 1 5 t 1 iP) , 



(22) 
(23) 



n 



(M - m ) 2 + p 2 
AG 




+ Tln 1 + e 




(24) 



where the energy Et- is defined by 



E± = ^{E±y + p> 



(25) 



§ In fact, every power of the quantum fluctuating fields that arises from expanding the functional 
determinant gives an additional factor of 1/ V N c and so Eq. (1 1 9 I t is a convenient way of organizing a 
l/iV c -expansion 
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E ± = E±Sn, (26) 



E = vV + M 2 . (27) 
In the limit T — > 0, the thermodynamic potential (|2~i|) reduces to 

o = (M -5 + " 2 - W / + (, - w - ^> 



4G L 7 (2tt) 2 

+ + , (28) 

which is in agreement with the result of Ref . [35] • The values of M and p are found 
by minimizing the thermodynamic potential f2, that is by solving the following gap 
equations 

aM = °< < 29 > 

-=0. (30) 

Differentiating the effective potential (JM]) with respect to M and /?, we obtain the gap 
equations 

- ^ " w < M j { W [1 " n(E ' " "> " " (E > + p)1 

+ ^[l-n(-E p --M)-n(-B p - +/.)]} , (31) 

° = £ - 2JVrf / (^{f^ - " (£ ' + - " - + " )] 

+ ^r[l - HE," - (.) - n(B; + ».)]}, (32) 
where n is the Fermi-Dirac distribution, 

Taking the limit T — > 0, these equations reduce to those obtained by Ebert and 
Klimenko [38] : 

2G y (2tt) 3 \ EE+ EEp J 

P OAr f d 3 p (8(E+-n) , -m) 



The dispersion relations for the quasiparticles are determined by the zeros of the 
functional determinant in Eq. (|19p . One finds [38] 



E u = E~ -n, E d = E+ - M (36) 

E u =E; +t i, Eg = E++p. (37) 

It is easy to show that the dispersion relations for the u and d-quarks are always 
gapped, while the dispersion relations for the u and d-quarks can be gapped or 
ungapped depending on the values of p. In the chiral limit, it follows directly from 
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the gap equations (|3ip and (|32p that there are no nonzero values for p, Sp, and T 
for which M and p are nonzero simultaneously. In the pion-condensed phase the 
dispersion relations for u and <i-quarks then reduce to 

E u = ^{ P -Sp) 2 +p 2 - H (38) 
E d = ^( P + Sp) 2 +p 2 - M ■ (39) 
For the u-quark, the dispersion relation is gapped or ungapped according to 

p > p , gapped spectrum , (40) 

p = p , ungapped quadratic spectrum , (41) 

p < p , ungapped linear spectrum . (42) 

The possibility of an ungapped quadratic quark spectrum in the context of dense 
baryonic matter was first discussed in Ref. 63J. The gaplessness of the <i-quark is 
determined by the line that is defined implicitly by the equation p — y/ (Sp) 2 + p 2 , 
which follows from Eq. (|39p . The window in the pion-condensed phase where the 
fermionic excitations are gapless is coined gapless pion condensation [38j . The 
situation here is analogous to what happens in colour superconductivity, except that 
in this case it is Sp that dictates the onset of gaplessness: Sp — A, where A is the 
superconducting gap [6] [64]. Moreover, in colour superconductivity Sp is called a 
stress parameter because it gives rise to a mismatch between the Fermi surfaces of the 
u and the <i-quarks, which imposes an extra energy cost (stress) on the formation of 
Cooper pairs. As long as the stress parameter is small enough compared to A, BCS 
pairing can occur [64j . In the context of pion condensation, p plays the role as a stress 
parameter since the mismatch is between the Fermi surfaces of the u and <i-quarks. 
Similarly, if the stress parameter p is small enough relative to p, pion condensation 
can occur. We will return to this point below. 

The integrals appearing in Eqs. (|24|) . (|3T|) . and ([32]) are ultraviolet divergent and 
one may impose a three-dimensional UV cutoff A to regulate them. Alternatively, one 
can introduce a form factor [65j , which falls off for large momenta. This ensures that 
the integrals converge in the ultraviolet. Given an ultraviolet cutoff A, the coupling 
constant G and the quark mass mo, we can use the Dyson equation for the quark 
propagator to determine the value for the chiral condensate in the vacuum. The 
Dyson equation reads [53] : 

M = „ 10 + 4 W /|^=, (43) 

where Nf is the number of flavours and N c is the number of colours. In the remainder 
of the paper, we also set Nf — 2, and N c = 3. This is simply the gap equation ([3T]) in 
the vacuum, i.e. for p = Sp = 0. 



3. Phase diagram 

In this section, we study the phase diagram at zero and finite temperature without 
taking charge neutrality into account. Since we are not considering colour super- 
conductivity, we restrict ourselves to quark chemical potentials p < 350 MeV. Since 
we are using an ultraviolet cutoff of approximately 650 MeV (see below), one should 
not trust results for the chemical potentials p and Sp above approximately 400 MeV. 
We therefore do not consider values above p = 350 MeV and Sp — 400 MeV in the 
calculations. The equilibrium values of M and p are obtained by solving numerically 
the gap equations ([3T]) and ([32]) . 



Pion condensation in a two-flavour NJL model 



8 



3.1. Chiral Limit 

In the chiral limit, the current quark mass vanishes, m = 0. For the numerical 
calculations, we choose an ultraviolet cutoff A = 650.9 MeV and a coupling constant 
G = 5.04 (GeV)~ 2 . Solving the Dyson equation (|43j) . this gives a constituent quark 
mass in the vacuum of M = 309.9 MeV. 

We mentioned in the previous section that there are no values of p, 5fi, and T 
such that M and p are simultaneously nonzero. In other words, the possible solutions 
to the gap equations are a) M = p = 0, b) M ^ 0, p = and c) p ^ 0, M = 0. In case 
a), the full symmetry of the Lagrangian is intact, while in case b) the U(1)l x U(1)r 
symmetry is broken spontaneously down to U(1)l+r by the chiral condensate. The 
breaking of the [/(l)-symmetry gives rise to a conventional Goldstone mod^jj In case 
c) the pion condensate breaks parity as well the U(1)l x U(1)r symmetry down to 
U(1)aj 3 . The latter transforms the left-handed and right-handed flavour doublet as 
ipL — > e laTa ipL, ipR — > e~ laTa ipR- The breaking of the £/(l)-symmetry gives rise to a 
conventional Goldstone boson with a linear dispersion relation for small values of the 
three-momentum. The ground state is therefore a pion superfluid 0. The dispersion 
relations for the composite fields a and ir in Eqs. (fT6|) - (fT7|) are found by first expanding 
the effective action (|T9|) to second order in the quantum fluctuating field a and 7r j . This 
gives rise to a 4 x 4 fluctuation matrix T(w,p), where the solutions to det T(u>,p) = 
determine the dispersion relations u>(j>). Details can be found in Ref. [39 . At finite 
chemical potential there is a mixing between the fields ffi and ii2 and one of the linear 
combinations is massless |38) . 

The pion condensate of the two-flavour NJL model at zero temperature is shown in 
Fig- El It is worth noting that along the /x-axis, i.e for dp = 0, the effective potential no 
longer depends on M and p separately, but rather on the combination M 2 + p 2 |39|. 
The effective potential then has the usual mexican-hat shape with infinitely many 
equivalent vacua, where we can choose anyone we wish. The chiral condensate can be 
rotated into pseudoscalar condensates via the axial flavour transformations 

^_» e <«»7»T.^ ) a = 1,2, 3. (44) 

This is analogous to what happens in colour superconductivity, where diquark 
condensates are rotated into pseudoscalar diquark ones via rotations similar to those 
in Eq. (|4"4"|) [T2"] . At finite quark mass or finite isospin chemical potential, the system 
becomes unstable against developing a nonzero pion condensate. At vanishing isospin 
chemical potential, parity is conserved in QCD and to enforce this we choose p = 
along the line 5p = 0. Consequently, the chiral condensate is nonvanishing and chiral 
symmetry is broken along the line 5fx = 0. Fig. [T] shows that the charged pions 
condense for any nonzero value of the isospin chemical potential This result is in 
accordance with that of Ebert and Klimenko [38 and calculations in the linear sigma 
model at finite isospin chemical potential ^33, 34^. Finally, we notice that the phase 
transition from the pion condensed phase to the chirally symmetric phase is first order. 
Notice that the critical isospin chemical potential 5p, is decreasing as a function of p. 

In Fig. [21 we show the pion condensate as a function of quark chemical potential 
and isospin chemical potential at T — 150 MeV. The region of pion condensation is 

|j For fi = 8fi = 0, i.e. in the vacuum, there is a broken SC/(2)-symmetry which gives rise to the 
three massless pions in the usual manner. 

If At finite density, Lorentz invariance is broken and the number of broken generators need not be 
the same as the number of Goldstone bosons 66, 67 68 . Goldstone bosons with quadratic dispersion 
relations appear and the system is not a superfluid in this case. 
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Figure 1. Pion condensate in the chiral limit as a function of quark chemical 
potential /i and Sfi at zero temperature. 



smaller than at T = and the transition to a chirally symmetric phase is now second 
order everywhere. 
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Figure 2. Pion condensate in the chiral limit as a function of /i and 5fi for 
T = 150 MeV. 

In Fig.[3l we show the pion condensate in the chiral limit as a function of /x and T 
for fixed value of isospin chemical potential, 8fi = 200 MeV. As the temperature and 
the quark chemical potential increase, the region of pion condensation decreases. For 
T = 0, the transition is first order. There is a line of first-order transitions starting at 
T = which ends at a critical point given by T = 0.11 GeV and fi = 0.22 GeV. The 
transition is of second order for larger values of T and smaller values of fi. 

As mentioned above, the quark chemical potential induces a stress in the system 
and for sufficiently large values of \i it is no longer energetically favourable to form a 
Bose condensate of charged pions. This is clearly seen in Fig.[U where we have plotted 
the thermodynamic potential f2(p) minus fl(p — 0) for three different values of /i with 
Sfi = 200 MeV and T = 0. From this figure it is evident that the transition is first 
order. 
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Figure 3. Pion condensate in the chiral limit as a function quark chemical p 
potential and temperature T for Sfi = 200 MeV. 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



P (GeV) 

Figure 4. Sl(p) — Q(p = 0) in the chiral limit at T = as a function of p for 
5p = 200 MeV and three different values of fi. Solid line: 100 MeV, dashed line: 
275.87MeV, and dotted line: 350 MeV. 

3. 2. Physical point 

At the physical point, we choose the parameters to = 5.5 MeV, A = 650.9 MeV, and 
G = 5.04 (GeV) -2 . Solving the Dyson equation |43|) . this gives a constituent quark 
vacuum mass M — 325.2 MeV, and the model reproduces the pion mass of = 140 
MeV. 

Due to the nonzero current quark mass, the chiral condensate M will always be 
nonzero and so chiral symmetry is never restored. For sufficiently high temperature 
or chemical potentials, the chiral condensate goes towards zero since the temperature- 
independent value of mo becomes irrelevant. 

The two possible solutions of the gap equations are therefore a) p — and b) 
p / 0. In the case a) the full symmetry of the Lagrangian is intact, while in case b) 
parity and the U(1)l+r symmetry are spontaneously broken. In the latter case, there 
is a conventional Goldstone mode and the ground state is a pion superfluid. 

In Fig. [5l we show the pion condensate as a function of T and S/i for fi = 0. 
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Figure 5. Pion condensate for fi = as a function of T and Sfi. 

The transition is second order everywhere with mean-held critical exponents, which 
is in agreement with the analysis using chiral perturbation theory [551 HZ]- However, 
lattice calculations [53] suggest that the transition is first order for fi] large enough. 
The line of hrst-order transitions ends at a tricritical point and the line of second-order 
transitions extends down to T = 0. This discrepancy between mean- field calculations 
and lattice calculations is most likely due to shortcomings of mean-field theory itself. 
Lattice simulations also suggest that the transition to a Bosc-Einstein condensed state 
coincides with the deconfinement phase transition to a quark-gluon plasma. 

In Fig. (6[ we show the pion condensate as a function of \i and b\x at the physical 
point for T = 0. Along the axis fi — 0, the transition to the pion condensed phase 
occurs at Sfi c = = 70 MeV. The phase transition for /x = is second order. 

The transition remains second order for 5[A C smaller than approximately 80 MeV. For 
larger values of Sfi it turns into a first-order transition. This is also shown in Fig. [71 
where the solid curve shows a first-order transition ending at a critical point, while the 
dashed line indicates a second order transition. This figure is very similar to Fig. 3 of 
Ref. [43] . 




Figure 6. Pion condensate at the physical point as a function of quark chemical 
potential fi and Sfi at zero temperature. 
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The results for the pion condensate are in qualitative agreement with that 
obtained by Barducci et al |37J. Quantitative difference are due to different quark- 
antiquark interaction terms in the Lagrangian as well as different ways of regulating 
the loop integrals. 
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Figure 7. Phase diagram at the physical point at T = 0. The solid curve 
indicates a first-order transition which ends at a critical point. The dashed curve 
indicates a second-order transition. 

In Fig. [51 we show the chiral and pion condensates as functions of the isospin 
chemical potential for /i = T = 0. The chiral condensate is constant from 8fi = to 
5/j, « TOtt/2, after which it drops rapidly towards zero at large <5/i. This behaviour is in 
agreement with the lattice simulations of Kogut and Sinclair [23] and the sigma-model 
calculations of He et al [33 . We notice the onset of pion condensation at Sfi « 
and that the condensate increases rapidly thereafter. This illustrates the competition 
between the two condensates. 
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Figure 8. Pion condensate (dashed curve) and chiral condensate (solid curve) of 
the two-flavour NJL model at the physical point as a function of 8fi at /i = T = 0. 



In Fig. [9l we show the pion condensate as a function of /i and S/i at the 
physical point for T = 150 MeV. The region of Bose condensation becomes smaller 
as the temperature increases, as expected. The phase transition is now second order 
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Figure 9. Pion condensate of the two-flavour NJL model at the physical point 
as a function of fi and Sfi for T = 150MeV. 



everywhere. 

In Fig. [TO] we show the phase diagram of the two-flavour NJL model as a function 
of T and 5/j, for (i = 0. The solid line is the chiral limit and the dashed line is the at 
the physical point. We note that the two curves are approaching each other for large 
values of 5ji, as expected. 
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Figure 10. Phase diagram as a function of Sfi and T for fi = 0. The solid line 
is the chiral limit and the dashed line is at the physical point. 

In Fig. 111! we show the chiral condensate as a function of fi and Sfi at the physical 
point and for T = 0. The chiral condensate goes to zero as the chemical potentials 
become large. 

4. Electric Charge Neutrality and /3-equilibrium 

In the previous section, we have calculated the phase diagram by finding the solutions 
to the gap equations and (|3"0|) . Dense matter inside stars should be neutral with 
respect to electric as well as colour charge, otherwise one would pay an enormous 
energy price [571 HE]- Since we are not considering colour superconducting phases, 
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Figure 11. Chiral condensate as a function of quark chemical potential fj, and 
Sfi for zero temperature. 



colour neutrality is automatically satisfied, while we have to impose electric charge 
neutrality. 

In addition to charge neutrality, matter should also be in (3 equilibrium, that is, 
weak-interaction processes such as 

u «-> d + e + + v , (45) 

should go with the same rate in both directions. If we assume that the neutrinos 
can leave the system, their chemical potential ll v vanishes. In chemical equlibrium, 
Eq. (|4"5|) then implies 

Mu = Md - Me • (46) 

The quark chemical potentials ll u and Lid, and the electron chemical potential ii e 
can be written in terms of the quark chemical potential fi and the electric chemical 
potential liq as 

2 

M« = M + gMQ i (47) 

Md = M - ^MQ , (48) 
Me = - MQ . (49) 

and so the system can be described in terms of the two independent chemical potentials 
fx and liq. In order to impose the constraint of charge neutrality, we additionally 
require that 

The constraint (|50|) implies that there is only one independent chemical potential, for 
example [i. Once we have picked a value for fx, the solutions to the gap equations ([29f 
and (|30|) and the neutrality constraint (|50|) determine the chiral condensate M . the 
charged pion condensate p, and the electric chemical potential liq. 

Note that the chemical potential appearing in Eq. (|24| is half the the sum of the 
quark chemical potentials ji u and Lid, and so according to Eqs. (|47[) and (|48[) . we need 
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to make the substitution p — > p, = p + Pq/6. In the remainder, we also replace Sp by 
Pq/2, which follows from Eqs. ([I5j). (|^7| . and (ggj) . 

In the following, we describe the electrons by a noninteracting Fermi gas. We 
then add to the Lagrangian (O, the term 

^electrons = 4>e (l^p + 7°^ ~ m e ) 1p e , (51) 

where ij) e denotes the electron field, e is the electron charge, and m e is the mass of the 
electron. The thermodynamic potential for the electrons is 

J3„ 



n, 



lectr 



- 2 

-Tin 



1 +e -/3(£ P +MQ) | 



1 _|_ e -P(Ep-HQ) 



(52) 



where £J p = \Jp 2 + mf. In the case of massless electrons, one can evaluate the integrals 
in Eq. (|5"2"|) exactly and one finds: 



n, 



4 ^Q T2 7tt 2 



lectrons- ^ g ^ . 

In the remainder, we neglect the electron mass. Adding the electron contribution 
Eq (|53p to Eq. ([24]) and differentiating with respect to pq, we obtain 







1 



^2 + q^Q T 
E+ 



1 - 3 



E~ 



d 3 P 
(2tt) 3 
1 



3— -3- 



£ 



-1 + 3 



1 



e 0(E++p.) + l 
1 



e P{E+-p.) + l 

1 + 3— I fl 3— , 

In the zero-temperature limit, Eq. (|54[) reduces to that obtained by Ebert and 
Klimenko I5UI: 



(54) 



o = -i 

3tt 2 



d 3 p 
(2^)3 



+ 3fl(£ p + - ji)^ - 39{E~ - /I) |=7 



T 



(55) 



5. Phase diagram Revisited 

In this section, we calculate the phase diagram of the two-flavour NJL model as 
a function of p and T imposing the electric charge neutrality constraint (|50[) . This 
equation and the gap equations (|3ip and (|32[) arc then solved simultaneously to obtain 
the equilibrium values of M and p for the neutral system. We are using the same 
parameter values as in the previous section. 



5.1. Chiral Limit 

Again the solutions to the gap equations are a) M = p = 0, b) M ^ 0, p = and c) 
p ± 0, M = 0. 

In Fig. 1121 we show the pion condensate as a function of quark chemical potential 
p and temperature T for neutral matter. The results at T — agree with those of 
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Ebert and Klimenko [3S]. Notice the black wedge starting in the corner \i = T = 0. 
In this area the electric chemical potential /j,q vanishes, giving rise to a nonzero chiral 
condensat43- For T = 0, the chiral condensate vanishes for quark chemical potentials 
larger that the a critical value of /Ui c = 297 MeV. For a quark chemical potential 
satisfying p,\ c < p, < /i2c, where ^2c = 329 MeV, the pion condensate is nonvanishing. 
For chemical potentials larger than p2c, the system is in the normal phase, where both 
condensates vanish. The robustness of this result is discussed in Sec. [6] Finally, we 
we notice that the transition from the pion-condensed phase to the symmetric phase 
is first order for T = 0. This first-order line ends at a critical point and the transition 
is second order all the way to p = 0. 



mo 
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Figure 12. Pion condensate in the chiral limit for neutral matter as a function 
of quark chemical potential fj, and temperature T. 



5. 2. Physical point 

Again the two possible solutions of the gap equations are 1) p = and 2) p ^ 0. It 
turns out that the only solution is p = 0, i.e. there is no charged pion condensate at 
the physical point. In other words, the isospin chemical potential pi — fiQ is always 
smaller than the pion mass. 

This is in accordance with the results of a recent study by Abuki et al [13] , where 
they investigated pion condensation in neutral matter as a function of the pion mass 
from the chiral limit all the way to the physical value of 139 MeV. They found a tiny 
window of pion condensation for pion masses below approximately 10 KeV. Thus pion 
condensation is very sensitive to the explicit chiral symmetry breaking due to a finite 
quark mass Too- 

In Fig. 1131 we show the chiral condensate for neutral matter as a function of 
temperature and quark chemical potential. The chiral condensate decreases with 
increasing fi and T, but never vanishes. 

+ Again the presence of the chiral condensate in the region where fiQ = is due to the fact that the 
effective potential depends on M 2 + p 2 (mexican-hat type potential) and we choose M ^ so that 
parity is unbroken in the vacuum. 
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Figure 13. Chiral condensate at the physical point for neutral matter as a 
function of temperature T and quark chemical potential fi 

6. Summary 

In this paper, we have calculated the phase diagram of the two-flavour NJL model 
as a function of the quark and the isospin chemical potential, and the temperature 
in the chiral limit and at the physical point. In the chiral limit, we have reproduced 
the results in Refs. [3H1I3S] at T = and generalized them to finite temperature. At 
the physical point, we get similar results for the phase structure as those obtained 
in Ref. [37]. The qualitative differences are due to different interaction terms in the 
Lagrangian and a different way of regulating the loop diagrams in the gap equations. 

It is natural to ask what happens at larger values of isospin chemical potential. We 
have extended our calculations of Sec. [3] for fi — ttiq — up to Sfi ~ 900 MeV. These 
calculations seem to indicate that there is a phase transition from a pion condensed 
phase to a chirally symmetric phase. For zero temperature, this transition is located 
at Sfi ~ 880 MeV. However, this result should not be trusted since the UV cutoff is 
651 MeV. In fact, there are reasons to believe that there is no phase transition as one 
increases Sfi. On general grounds, one can show that there is a nonzero condensate 
of the form (J75 it) for larger values of hi [55]. Thus one expects a BEC-BCS type of 
crossover as the quarks bound in the pions become weakly bound due to the fact that 
the QCD running coupling becomes weaker with increasing chemical potential /ij. 

Our main result comes from imposing the constraints of electric charge neutrality 
and weak equilibrium. In Ref. [39], the authors calculate the phases of neutral matter 
with another set of parameters at zero temperature. Using an ultraviolet cutoff of 600 
MeV, a coupling constant G = 6.82 (GeV) -2 , and a constituent quark mass of 400 
MeV, their numerical analysis shows that the phase structure differs from the first 
set of parameters. At T — 0, there is in this case no phase with a pion condensate, 
but a phase transition directly from a phase of broken chiral symmetry to a chirally 
symmetric phase at a critical quark chemical potential of ji c = 386.2 MeV. This seems 
to indicate that the window for pion condensation at T — that can be seen in Fig. [T^] 
is not a robust result. 

We have seen that the charge neutrality constraint changes the phase diagram 
since it rules out charged pion condensation at the physical point. This is true for pion 
masses larger than approximately 10 KeV [43]. However, in the presence of a dense 
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neutrino gas, pion condensation again becomes a possibility even with realistic values 
of m v , as shown in Ref. [IB]. In nature this situation arises in supernova explosions 
and possibly at the early stages of the evolution of a neutron star. It does not arise 
in stable matter such as we have considered here, because the weakly interacting 
neutrinos will have had time to leave the system. 

A natural next step would be to investigate the competition between pion 
and kaon condensation in neutral matter using a three-flavour NJL model at finite 
temperature and finite quark chemical potentials (j, u , fid, and /i s . Calculations without 
the neutrality constraint have already been done by Barducci et al |69j . 
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